\(\int \frac {\cot ^3(d+e x)}{\sqrt {a+b \tan ^2(d+e x)+c \tan ^4(d+e x)}} \, dx\) [40]

   Optimal result
   Rubi [A] (verified)
   Mathematica [A] (verified)
   Maple [F]
   Fricas [A] (verification not implemented)
   Sympy [F]
   Maxima [F(-1)]
   Giac [F]
   Mupad [F(-1)]

Optimal result

Integrand size = 35, antiderivative size = 249 \[ \int \frac {\cot ^3(d+e x)}{\sqrt {a+b \tan ^2(d+e x)+c \tan ^4(d+e x)}} \, dx=\frac {\text {arctanh}\left (\frac {2 a+b \tan ^2(d+e x)}{2 \sqrt {a} \sqrt {a+b \tan ^2(d+e x)+c \tan ^4(d+e x)}}\right )}{2 \sqrt {a} e}+\frac {b \text {arctanh}\left (\frac {2 a+b \tan ^2(d+e x)}{2 \sqrt {a} \sqrt {a+b \tan ^2(d+e x)+c \tan ^4(d+e x)}}\right )}{4 a^{3/2} e}-\frac {\text {arctanh}\left (\frac {2 a-b+(b-2 c) \tan ^2(d+e x)}{2 \sqrt {a-b+c} \sqrt {a+b \tan ^2(d+e x)+c \tan ^4(d+e x)}}\right )}{2 \sqrt {a-b+c} e}-\frac {\cot ^2(d+e x) \sqrt {a+b \tan ^2(d+e x)+c \tan ^4(d+e x)}}{2 a e} \]

[Out]

1/4*b*arctanh(1/2*(2*a+b*tan(e*x+d)^2)/a^(1/2)/(a+b*tan(e*x+d)^2+c*tan(e*x+d)^4)^(1/2))/a^(3/2)/e+1/2*arctanh(
1/2*(2*a+b*tan(e*x+d)^2)/a^(1/2)/(a+b*tan(e*x+d)^2+c*tan(e*x+d)^4)^(1/2))/e/a^(1/2)-1/2*arctanh(1/2*(2*a-b+(b-
2*c)*tan(e*x+d)^2)/(a-b+c)^(1/2)/(a+b*tan(e*x+d)^2+c*tan(e*x+d)^4)^(1/2))/e/(a-b+c)^(1/2)-1/2*cot(e*x+d)^2*(a+
b*tan(e*x+d)^2+c*tan(e*x+d)^4)^(1/2)/a/e

Rubi [A] (verified)

Time = 0.37 (sec) , antiderivative size = 249, normalized size of antiderivative = 1.00, number of steps used = 11, number of rules used = 6, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.171, Rules used = {3781, 1265, 974, 744, 738, 212} \[ \int \frac {\cot ^3(d+e x)}{\sqrt {a+b \tan ^2(d+e x)+c \tan ^4(d+e x)}} \, dx=\frac {b \text {arctanh}\left (\frac {2 a+b \tan ^2(d+e x)}{2 \sqrt {a} \sqrt {a+b \tan ^2(d+e x)+c \tan ^4(d+e x)}}\right )}{4 a^{3/2} e}+\frac {\text {arctanh}\left (\frac {2 a+b \tan ^2(d+e x)}{2 \sqrt {a} \sqrt {a+b \tan ^2(d+e x)+c \tan ^4(d+e x)}}\right )}{2 \sqrt {a} e}-\frac {\text {arctanh}\left (\frac {2 a+(b-2 c) \tan ^2(d+e x)-b}{2 \sqrt {a-b+c} \sqrt {a+b \tan ^2(d+e x)+c \tan ^4(d+e x)}}\right )}{2 e \sqrt {a-b+c}}-\frac {\cot ^2(d+e x) \sqrt {a+b \tan ^2(d+e x)+c \tan ^4(d+e x)}}{2 a e} \]

[In]

Int[Cot[d + e*x]^3/Sqrt[a + b*Tan[d + e*x]^2 + c*Tan[d + e*x]^4],x]

[Out]

ArcTanh[(2*a + b*Tan[d + e*x]^2)/(2*Sqrt[a]*Sqrt[a + b*Tan[d + e*x]^2 + c*Tan[d + e*x]^4])]/(2*Sqrt[a]*e) + (b
*ArcTanh[(2*a + b*Tan[d + e*x]^2)/(2*Sqrt[a]*Sqrt[a + b*Tan[d + e*x]^2 + c*Tan[d + e*x]^4])])/(4*a^(3/2)*e) -
ArcTanh[(2*a - b + (b - 2*c)*Tan[d + e*x]^2)/(2*Sqrt[a - b + c]*Sqrt[a + b*Tan[d + e*x]^2 + c*Tan[d + e*x]^4])
]/(2*Sqrt[a - b + c]*e) - (Cot[d + e*x]^2*Sqrt[a + b*Tan[d + e*x]^2 + c*Tan[d + e*x]^4])/(2*a*e)

Rule 212

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1/(Rt[a, 2]*Rt[-b, 2]))*ArcTanh[Rt[-b, 2]*(x/Rt[a, 2])], x]
 /; FreeQ[{a, b}, x] && NegQ[a/b] && (GtQ[a, 0] || LtQ[b, 0])

Rule 738

Int[1/(((d_.) + (e_.)*(x_))*Sqrt[(a_.) + (b_.)*(x_) + (c_.)*(x_)^2]), x_Symbol] :> Dist[-2, Subst[Int[1/(4*c*d
^2 - 4*b*d*e + 4*a*e^2 - x^2), x], x, (2*a*e - b*d - (2*c*d - b*e)*x)/Sqrt[a + b*x + c*x^2]], x] /; FreeQ[{a,
b, c, d, e}, x] && NeQ[b^2 - 4*a*c, 0] && NeQ[2*c*d - b*e, 0]

Rule 744

Int[((d_.) + (e_.)*(x_))^(m_)*((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)^(p_), x_Symbol] :> Simp[e*(d + e*x)^(m + 1)*
((a + b*x + c*x^2)^(p + 1)/((m + 1)*(c*d^2 - b*d*e + a*e^2))), x] + Dist[(2*c*d - b*e)/(2*(c*d^2 - b*d*e + a*e
^2)), Int[(d + e*x)^(m + 1)*(a + b*x + c*x^2)^p, x], x] /; FreeQ[{a, b, c, d, e, m, p}, x] && NeQ[b^2 - 4*a*c,
 0] && NeQ[c*d^2 - b*d*e + a*e^2, 0] && NeQ[2*c*d - b*e, 0] && EqQ[m + 2*p + 3, 0]

Rule 974

Int[((d_.) + (e_.)*(x_))^(m_)*((f_.) + (g_.)*(x_))^(n_)*((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)^(p_.), x_Symbol] :
> Int[ExpandIntegrand[(d + e*x)^m*(f + g*x)^n*(a + b*x + c*x^2)^p, x], x] /; FreeQ[{a, b, c, d, e, f, g}, x] &
& NeQ[e*f - d*g, 0] && NeQ[b^2 - 4*a*c, 0] && NeQ[c*d^2 - b*d*e + a*e^2, 0] && (IntegerQ[p] || (ILtQ[m, 0] &&
ILtQ[n, 0])) &&  !(IGtQ[m, 0] || IGtQ[n, 0])

Rule 1265

Int[(x_)^(m_.)*((d_) + (e_.)*(x_)^2)^(q_.)*((a_) + (b_.)*(x_)^2 + (c_.)*(x_)^4)^(p_.), x_Symbol] :> Dist[1/2,
Subst[Int[x^((m - 1)/2)*(d + e*x)^q*(a + b*x + c*x^2)^p, x], x, x^2], x] /; FreeQ[{a, b, c, d, e, p, q}, x] &&
 IntegerQ[(m - 1)/2]

Rule 3781

Int[tan[(d_.) + (e_.)*(x_)]^(m_.)*((a_.) + (b_.)*((f_.)*tan[(d_.) + (e_.)*(x_)])^(n_.) + (c_.)*((f_.)*tan[(d_.
) + (e_.)*(x_)])^(n2_.))^(p_), x_Symbol] :> Dist[f/e, Subst[Int[(x/f)^m*((a + b*x^n + c*x^(2*n))^p/(f^2 + x^2)
), x], x, f*Tan[d + e*x]], x] /; FreeQ[{a, b, c, d, e, f, m, n, p}, x] && EqQ[n2, 2*n] && NeQ[b^2 - 4*a*c, 0]

Rubi steps \begin{align*} \text {integral}& = \frac {\text {Subst}\left (\int \frac {1}{x^3 \left (1+x^2\right ) \sqrt {a+b x^2+c x^4}} \, dx,x,\tan (d+e x)\right )}{e} \\ & = \frac {\text {Subst}\left (\int \frac {1}{x^2 (1+x) \sqrt {a+b x+c x^2}} \, dx,x,\tan ^2(d+e x)\right )}{2 e} \\ & = \frac {\text {Subst}\left (\int \left (\frac {1}{x^2 \sqrt {a+b x+c x^2}}-\frac {1}{x \sqrt {a+b x+c x^2}}+\frac {1}{(1+x) \sqrt {a+b x+c x^2}}\right ) \, dx,x,\tan ^2(d+e x)\right )}{2 e} \\ & = \frac {\text {Subst}\left (\int \frac {1}{x^2 \sqrt {a+b x+c x^2}} \, dx,x,\tan ^2(d+e x)\right )}{2 e}-\frac {\text {Subst}\left (\int \frac {1}{x \sqrt {a+b x+c x^2}} \, dx,x,\tan ^2(d+e x)\right )}{2 e}+\frac {\text {Subst}\left (\int \frac {1}{(1+x) \sqrt {a+b x+c x^2}} \, dx,x,\tan ^2(d+e x)\right )}{2 e} \\ & = -\frac {\cot ^2(d+e x) \sqrt {a+b \tan ^2(d+e x)+c \tan ^4(d+e x)}}{2 a e}+\frac {\text {Subst}\left (\int \frac {1}{4 a-x^2} \, dx,x,\frac {2 a+b \tan ^2(d+e x)}{\sqrt {a+b \tan ^2(d+e x)+c \tan ^4(d+e x)}}\right )}{e}-\frac {\text {Subst}\left (\int \frac {1}{4 a-4 b+4 c-x^2} \, dx,x,\frac {2 a-b-(-b+2 c) \tan ^2(d+e x)}{\sqrt {a+b \tan ^2(d+e x)+c \tan ^4(d+e x)}}\right )}{e}-\frac {b \text {Subst}\left (\int \frac {1}{x \sqrt {a+b x+c x^2}} \, dx,x,\tan ^2(d+e x)\right )}{4 a e} \\ & = \frac {\text {arctanh}\left (\frac {2 a+b \tan ^2(d+e x)}{2 \sqrt {a} \sqrt {a+b \tan ^2(d+e x)+c \tan ^4(d+e x)}}\right )}{2 \sqrt {a} e}-\frac {\text {arctanh}\left (\frac {2 a-b+(b-2 c) \tan ^2(d+e x)}{2 \sqrt {a-b+c} \sqrt {a+b \tan ^2(d+e x)+c \tan ^4(d+e x)}}\right )}{2 \sqrt {a-b+c} e}-\frac {\cot ^2(d+e x) \sqrt {a+b \tan ^2(d+e x)+c \tan ^4(d+e x)}}{2 a e}+\frac {b \text {Subst}\left (\int \frac {1}{4 a-x^2} \, dx,x,\frac {2 a+b \tan ^2(d+e x)}{\sqrt {a+b \tan ^2(d+e x)+c \tan ^4(d+e x)}}\right )}{2 a e} \\ & = \frac {\text {arctanh}\left (\frac {2 a+b \tan ^2(d+e x)}{2 \sqrt {a} \sqrt {a+b \tan ^2(d+e x)+c \tan ^4(d+e x)}}\right )}{2 \sqrt {a} e}+\frac {b \text {arctanh}\left (\frac {2 a+b \tan ^2(d+e x)}{2 \sqrt {a} \sqrt {a+b \tan ^2(d+e x)+c \tan ^4(d+e x)}}\right )}{4 a^{3/2} e}-\frac {\text {arctanh}\left (\frac {2 a-b+(b-2 c) \tan ^2(d+e x)}{2 \sqrt {a-b+c} \sqrt {a+b \tan ^2(d+e x)+c \tan ^4(d+e x)}}\right )}{2 \sqrt {a-b+c} e}-\frac {\cot ^2(d+e x) \sqrt {a+b \tan ^2(d+e x)+c \tan ^4(d+e x)}}{2 a e} \\ \end{align*}

Mathematica [A] (verified)

Time = 4.42 (sec) , antiderivative size = 188, normalized size of antiderivative = 0.76 \[ \int \frac {\cot ^3(d+e x)}{\sqrt {a+b \tan ^2(d+e x)+c \tan ^4(d+e x)}} \, dx=\frac {(2 a+b) \text {arctanh}\left (\frac {2 a+b \tan ^2(d+e x)}{2 \sqrt {a} \sqrt {a+b \tan ^2(d+e x)+c \tan ^4(d+e x)}}\right )+2 \sqrt {a} \left (-\frac {a \text {arctanh}\left (\frac {2 a-b+(b-2 c) \tan ^2(d+e x)}{2 \sqrt {a-b+c} \sqrt {a+b \tan ^2(d+e x)+c \tan ^4(d+e x)}}\right )}{\sqrt {a-b+c}}-\cot ^2(d+e x) \sqrt {a+b \tan ^2(d+e x)+c \tan ^4(d+e x)}\right )}{4 a^{3/2} e} \]

[In]

Integrate[Cot[d + e*x]^3/Sqrt[a + b*Tan[d + e*x]^2 + c*Tan[d + e*x]^4],x]

[Out]

((2*a + b)*ArcTanh[(2*a + b*Tan[d + e*x]^2)/(2*Sqrt[a]*Sqrt[a + b*Tan[d + e*x]^2 + c*Tan[d + e*x]^4])] + 2*Sqr
t[a]*(-((a*ArcTanh[(2*a - b + (b - 2*c)*Tan[d + e*x]^2)/(2*Sqrt[a - b + c]*Sqrt[a + b*Tan[d + e*x]^2 + c*Tan[d
 + e*x]^4])])/Sqrt[a - b + c]) - Cot[d + e*x]^2*Sqrt[a + b*Tan[d + e*x]^2 + c*Tan[d + e*x]^4]))/(4*a^(3/2)*e)

Maple [F]

\[\int \frac {\cot \left (e x +d \right )^{3}}{\sqrt {a +b \tan \left (e x +d \right )^{2}+c \tan \left (e x +d \right )^{4}}}d x\]

[In]

int(cot(e*x+d)^3/(a+b*tan(e*x+d)^2+c*tan(e*x+d)^4)^(1/2),x)

[Out]

int(cot(e*x+d)^3/(a+b*tan(e*x+d)^2+c*tan(e*x+d)^4)^(1/2),x)

Fricas [A] (verification not implemented)

none

Time = 1.29 (sec) , antiderivative size = 1350, normalized size of antiderivative = 5.42 \[ \int \frac {\cot ^3(d+e x)}{\sqrt {a+b \tan ^2(d+e x)+c \tan ^4(d+e x)}} \, dx=\text {Too large to display} \]

[In]

integrate(cot(e*x+d)^3/(a+b*tan(e*x+d)^2+c*tan(e*x+d)^4)^(1/2),x, algorithm="fricas")

[Out]

[1/8*(2*sqrt(a - b + c)*a^2*log(((b^2 + 4*(a - 2*b)*c + 8*c^2)*tan(e*x + d)^4 + 2*(4*a*b - 3*b^2 - 4*(a - b)*c
)*tan(e*x + d)^2 - 4*sqrt(c*tan(e*x + d)^4 + b*tan(e*x + d)^2 + a)*((b - 2*c)*tan(e*x + d)^2 + 2*a - b)*sqrt(a
 - b + c) + 8*a^2 - 8*a*b + b^2 + 4*a*c)/(tan(e*x + d)^4 + 2*tan(e*x + d)^2 + 1))*tan(e*x + d)^2 + (2*a^2 - a*
b - b^2 + (2*a + b)*c)*sqrt(a)*log(((b^2 + 4*a*c)*tan(e*x + d)^4 + 8*a*b*tan(e*x + d)^2 + 4*sqrt(c*tan(e*x + d
)^4 + b*tan(e*x + d)^2 + a)*(b*tan(e*x + d)^2 + 2*a)*sqrt(a) + 8*a^2)/tan(e*x + d)^4)*tan(e*x + d)^2 - 4*sqrt(
c*tan(e*x + d)^4 + b*tan(e*x + d)^2 + a)*(a^2 - a*b + a*c))/((a^3 - a^2*b + a^2*c)*e*tan(e*x + d)^2), 1/4*(sqr
t(a - b + c)*a^2*log(((b^2 + 4*(a - 2*b)*c + 8*c^2)*tan(e*x + d)^4 + 2*(4*a*b - 3*b^2 - 4*(a - b)*c)*tan(e*x +
 d)^2 - 4*sqrt(c*tan(e*x + d)^4 + b*tan(e*x + d)^2 + a)*((b - 2*c)*tan(e*x + d)^2 + 2*a - b)*sqrt(a - b + c) +
 8*a^2 - 8*a*b + b^2 + 4*a*c)/(tan(e*x + d)^4 + 2*tan(e*x + d)^2 + 1))*tan(e*x + d)^2 - (2*a^2 - a*b - b^2 + (
2*a + b)*c)*sqrt(-a)*arctan(1/2*sqrt(c*tan(e*x + d)^4 + b*tan(e*x + d)^2 + a)*(b*tan(e*x + d)^2 + 2*a)*sqrt(-a
)/(a*c*tan(e*x + d)^4 + a*b*tan(e*x + d)^2 + a^2))*tan(e*x + d)^2 - 2*sqrt(c*tan(e*x + d)^4 + b*tan(e*x + d)^2
 + a)*(a^2 - a*b + a*c))/((a^3 - a^2*b + a^2*c)*e*tan(e*x + d)^2), -1/8*(4*a^2*sqrt(-a + b - c)*arctan(-1/2*sq
rt(c*tan(e*x + d)^4 + b*tan(e*x + d)^2 + a)*((b - 2*c)*tan(e*x + d)^2 + 2*a - b)*sqrt(-a + b - c)/(((a - b)*c
+ c^2)*tan(e*x + d)^4 + (a*b - b^2 + b*c)*tan(e*x + d)^2 + a^2 - a*b + a*c))*tan(e*x + d)^2 - (2*a^2 - a*b - b
^2 + (2*a + b)*c)*sqrt(a)*log(((b^2 + 4*a*c)*tan(e*x + d)^4 + 8*a*b*tan(e*x + d)^2 + 4*sqrt(c*tan(e*x + d)^4 +
 b*tan(e*x + d)^2 + a)*(b*tan(e*x + d)^2 + 2*a)*sqrt(a) + 8*a^2)/tan(e*x + d)^4)*tan(e*x + d)^2 + 4*sqrt(c*tan
(e*x + d)^4 + b*tan(e*x + d)^2 + a)*(a^2 - a*b + a*c))/((a^3 - a^2*b + a^2*c)*e*tan(e*x + d)^2), -1/4*(2*a^2*s
qrt(-a + b - c)*arctan(-1/2*sqrt(c*tan(e*x + d)^4 + b*tan(e*x + d)^2 + a)*((b - 2*c)*tan(e*x + d)^2 + 2*a - b)
*sqrt(-a + b - c)/(((a - b)*c + c^2)*tan(e*x + d)^4 + (a*b - b^2 + b*c)*tan(e*x + d)^2 + a^2 - a*b + a*c))*tan
(e*x + d)^2 + (2*a^2 - a*b - b^2 + (2*a + b)*c)*sqrt(-a)*arctan(1/2*sqrt(c*tan(e*x + d)^4 + b*tan(e*x + d)^2 +
 a)*(b*tan(e*x + d)^2 + 2*a)*sqrt(-a)/(a*c*tan(e*x + d)^4 + a*b*tan(e*x + d)^2 + a^2))*tan(e*x + d)^2 + 2*sqrt
(c*tan(e*x + d)^4 + b*tan(e*x + d)^2 + a)*(a^2 - a*b + a*c))/((a^3 - a^2*b + a^2*c)*e*tan(e*x + d)^2)]

Sympy [F]

\[ \int \frac {\cot ^3(d+e x)}{\sqrt {a+b \tan ^2(d+e x)+c \tan ^4(d+e x)}} \, dx=\int \frac {\cot ^{3}{\left (d + e x \right )}}{\sqrt {a + b \tan ^{2}{\left (d + e x \right )} + c \tan ^{4}{\left (d + e x \right )}}}\, dx \]

[In]

integrate(cot(e*x+d)**3/(a+b*tan(e*x+d)**2+c*tan(e*x+d)**4)**(1/2),x)

[Out]

Integral(cot(d + e*x)**3/sqrt(a + b*tan(d + e*x)**2 + c*tan(d + e*x)**4), x)

Maxima [F(-1)]

Timed out. \[ \int \frac {\cot ^3(d+e x)}{\sqrt {a+b \tan ^2(d+e x)+c \tan ^4(d+e x)}} \, dx=\text {Timed out} \]

[In]

integrate(cot(e*x+d)^3/(a+b*tan(e*x+d)^2+c*tan(e*x+d)^4)^(1/2),x, algorithm="maxima")

[Out]

Timed out

Giac [F]

\[ \int \frac {\cot ^3(d+e x)}{\sqrt {a+b \tan ^2(d+e x)+c \tan ^4(d+e x)}} \, dx=\int { \frac {\cot \left (e x + d\right )^{3}}{\sqrt {c \tan \left (e x + d\right )^{4} + b \tan \left (e x + d\right )^{2} + a}} \,d x } \]

[In]

integrate(cot(e*x+d)^3/(a+b*tan(e*x+d)^2+c*tan(e*x+d)^4)^(1/2),x, algorithm="giac")

[Out]

sage0*x

Mupad [F(-1)]

Timed out. \[ \int \frac {\cot ^3(d+e x)}{\sqrt {a+b \tan ^2(d+e x)+c \tan ^4(d+e x)}} \, dx=\int \frac {{\mathrm {cot}\left (d+e\,x\right )}^3}{\sqrt {c\,{\mathrm {tan}\left (d+e\,x\right )}^4+b\,{\mathrm {tan}\left (d+e\,x\right )}^2+a}} \,d x \]

[In]

int(cot(d + e*x)^3/(a + b*tan(d + e*x)^2 + c*tan(d + e*x)^4)^(1/2),x)

[Out]

int(cot(d + e*x)^3/(a + b*tan(d + e*x)^2 + c*tan(d + e*x)^4)^(1/2), x)